Non-reciprocal optical reflection from a bidimensional array of 
subwavelength holes in a metallic film 
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Using simulations and theoretical arguments we investigate the specular reflection of a perforated 
gold film deposited on a glass substrate. A square lattice of cylindrical holes is assumed to produce 
the periodic lateral corrugation needed to hybridize the surface plasmons with radiative states. It 
is shown that, contrasting transmission approaches, a knowledge of the reflection on either side 
of the fllm provides separate information on the gold-vacuum surface plasmons and on the gold- 
glass interface plasmons. Recent experimental data on a specific implementation of this system are 
reexamined; these show a good agreement between the measured refiections and the simulations in 
both directions of incident wave probes. This confirms the importance of taking into account the 
reflection asymmetry in the far-field assessment of surface plasmons properties. 

PACS numbers: 78.20.-e, 42.79.Dj, 42.25.Bs, 42.25.Gy 



In recent years, subwavelength hole arrays in metal- 
lic films have received considerable attention, mainly be- 
cause of the experimental findings of Thomas Ebbeseiii. 
Such systems show a particularly intriguing zeroth order 
optical transmission : first, the specularly transmitted 
beam shows an intensity higher than that of the incident 
beam filtered by the geometrical aperture of the holes; 
second, the transmitted beam develops a strong charac- 
teristic wavelength dependencei. Recent theoretical and 
experimental studies have tried to clarify the detailed 
mechanisms involved in these phenomenaiSiSiiiSiSiLS. In 
this context, it was rapidly recognized that surface plas- 
mons (SPs) could be a key to the understanding of the 
transmission measurement3iiSi2i^iSiSiL&. 

To be more precise, when the incident light illuminates 
the array of apertures, many diffraction orders are gen- 
erated. Some of them are evanescent and exhibit a res- 
onant coupling with SPs. They contribute to the ze- 
roth diffraction order as a result of multiscattering and 
the strong exaltation of the electromagnetic field asso- 
ciated with resonant processes results in an enhanced 
transmissioniiSiSi^iSiSiLS. In a recent worls& we have sug- 
gested redefining the role of SPs in Ebbesen experiments 
and placing it in the required context of resonant Wood 
anomaliesS. In so doing, we have shown that the trans- 
mission spectrum could be better depicted as a series 
of Fano profiles. These recognizable lineshapes result 
from the interference of non resonant transfers with res- 
onant transfers which involve the film eigenmodes and 
evanescent diffraction orders. We were then able to draw 
attention to the fact that each transmission maximum- 
minimum pair can be viewed as a Fano profileifi. On 
quantitative grounds, the interplay of such an assymetric 
Fano profile actually dissociates the presence of a maxi- 
mum or minimum in the transmission spectrum from the 
expected location of an eigenmode resonance^. 

It is somewhat surprising that reflection has not been 
as widely studied as transmission. One of our early con- 



tributions to this field described the theoretical refiection 
of an array of holes in a chromium film^ and, at about 
the same time, Altewischer et al reported experimen- 
tal reflection studies which emphasized its non-reciprocal 
propertieai. As is well-known, there is a significant differ- 
ence between the reflection spectra for light beams com- 
ing from one side or the other side of the device. Such 
difference is not observed in the transmission spectra, 
in agreement with the reciprocity theoremiii. The dis- 
symmetry observed in the reflection spectra usually con- 
tains very signiflcant information and both reflection and 
transmission should be thoroughly understood before a 
full knowledge of the optical properties of the system can 
be claimedSii. The present paper proposes a theoretical 
investigation of the device studied by Altewischer et a&. 
This will provide an opportunity to establish the role of 
SPs in these experiments and, at the same time, relate 
the far-field transfer coefficients to the very local SPs dy- 
namics. 

The simulations in this paper are based on a coupled 
modes method which combines a scattering matrix for- 
malism with a plane wave representation of the fields-. 
The inherently slow convergence of this representation 
is here accelerated by a Li algorithm—. This technique 
provides a computation scheme for the amplitude and 
polarization (s or p) of reflected and transmitted fields 
in any diffracted ordei"^. 

FigCl describes the system under study, chosen 
to model the samples experimentally investigated by 
Altewischer—. These samples can be described as a square 
lattice grating of circular holes perforated in a thin gold 
film deposited on a glass substrate. In accordance with 
experiments, 100 nm is set for the radius of each hole 
and the lattice parameter is taken to be a = 700 nm. 
The film thickness is matched to the experimental value 
of 200 nm. The gold permittivity is obtained from ex- 
perimental tablesi^. The substrate slab is modelled by 
a semi-infinite glass medium with a refractive index of 
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FIG. 1: A view of the system under studj'^'. Transmission and 
reflection are calculated for zero diffraction order and normal 
incidence. 



FIG. 2: (color online) (a) : Calculated transmission (curve 
marked T) and reflection for zero diffraction order on the 
structure of Fig0 On reflection curves, the solid (dashed) 
line refers to incidence from the vacuum (glass) side. The 
reflection nonreciprocity (difference according to incident di- 
rections) takes the form of two peaks at wavelengths (1') 
and (2') for vacuum incidence, while only one at wavelength 
(2') remains in glass incidence. Other panels : (b) com- 
puted (dashed) transmission and, (c) and (d) reflections (from 
vacuum and glass, respectively) compared to measurements 
(solid) from Altewischer et al~. Note that experimental values 
are not available above 1100 nm. 

— 1.51—. Note that the experimental device exhib- 
ited a very thin bonding layer of titanium which separates 
the gold film from the substrate^. The vanishing thick- 
ness of this film justifies that it will not be considered 
in our calculations. The linearly polarized incident light 
is normal to the interface. Its electric-field vector is ori- 
ented along the [1, 0] direction, which joins the centers of 
nearest- neighboring holes. 

Fig|2Ia) shows the specular reflection as a function of 
the incident wavelength in both directions of propaga- 
tion, i.e. an incidence from vacuum or from glass. The 
transmission is also shown, with a 5-times expansion of 
the vertical scaled. The reflection presents strong vari- 
ations near wavelengths 757 nm (1'), 812 nm (2') and 
1146 nm (3'). These variations are fully consistent with 
experimental spectral features found at 747 nm, 810 nm 
and, presumably, at just over 1100 nm^. Fig[2Ib), |2Ic), 
and|2td) reproduces the curves shown on Figl^^a), but 
exposes them to a comparison with the corresponding ex- 
perimental data. Despite the necessary idealization step 
in designing the modelled theoretical device, we note an 
excellent agreement in all these comparisonai. In partic- 
ular, the relative locations of transmission maxima and 
reflection minima are correctly predicted. The only dis- 
crepancies which may still raise questions seem to be the 
weak oscillations found in all calculated spectra between 
850 nm and 1050 nm, not found in experiment. The spec- 
tral line shapes are reminiscent of those described for sim- 
ilar systems in recent computation^ and experimental^ 
reports. The difference between the reflections on both 
sides of the film, i.e. the reflection nonreciprocity, is evi- 
dent from the comparison of Figl^Jc) and Figl^^d) . These 
non-reciprocal properties of the reflection were also men- 
tioned by Barnes et al^ when describing their optical 
properties of structured metallic films. 

In the present investigation the coupling between the 
various evanescent p-polarized diffraction orders and the 
surface plasmons is of prime importance. Because of the 
periodic corrugation of the film, the diffracted surface 
plasmons are sampled by the diffracted probe field at 



FIG. 3: (color online) Electric fleld modulus spectrum for 
various diffraction orders for light impinging from the vac- 
uum side (solid line) or from the glass side (dashed line) of 
the film, (a) Diffracted field of orders (±1,0) on the sub- 
strate/metal interface, (b) Diffracted field of orders (±1,0) 
at the vacuum/metal interface, (c) Diffracted field of orders 
(±1,±1) at the substrate/metal interface. The amplitude of 
the incident field is taken to be 1 V.m~^. (d) Absorption from 
the incident wave against wavelength, for incidences from the 
vacuum side (solid line) and from the glass side (dashed line). 

many frequencies. The resonant coupling of surface plas- 
mons to p-polarized diffraction orders at matching fre- 
quencies and surface wave vectors will lead to what will 
be called "resonant diffraction orders". For an isolated 
surface plasmon which develops on a single interface, the 
wavelength which characterizes each resonant diffraction 
order is given hyi^ 

In this equation, X^p — 2ttc/uJsp where ujsp expresses the 
resonant angular frequency, a is the grating parameter 
and denotes the vector g of the reciprocal lattice 

such that g = (27r/a)(i, j). e^n is the metal permittivity 
and £„ is either the permittivity of glass (Sg) or vacuum 
(e. = 1). 

For the resonant orders (±1,0) we get at sub- 
strate/metal and vacuum/metal interfaces Asp™ ^^'^^ — 

1084.6 nm and A^p™ '''^'^^ — 720.5 nm respectively. Both 
wavelengths are indicated by vertical dashes (3) and 
(1) respectively in the Fig[21 and Fig|31 For the res- 
onant orders (±1,±1) at substrate/metal interface, we 
get Asp™ ^^'^^ = 787.7 nm indicated by dash (2). Dash 
(0) corresponds to the resonant orders (±2, 0) at sub- 
strate/metal interface for Asp™ = 608.1 nm. This 
resonance appears in transmission as the first peak at 
610 nm. 

The contrast between the refiection spectra for both 
incidence directions allows additional information to be 
obtained on processes involving specifically the surface 
plasmons at both interfacea^i. For light arriving from the 
vacuum, two dips at the wavelengths (1') and (2') are ob- 
served. By contrast, for light arriving from the glass side 
a significant dip shows up at the wavelength (3') while 
(1') almost disappears. Altewischer et a& suggested that 
such an asymmetry in the reflection spectra could lead 
to the identification of the interface on which the SPs are 
dominantly excited. We think that, furthermore, these 
results indicate that the excited resonances are mainly 
localized on the interface which is illuminated first and 
thus their effects are more easily perceived with refiection 
than with transmission. 

To confirm and detail this interpretation, we have com- 
puted the amplitudes of main resonant diffraction orders 
as a function of wavelength and illumination directions. 
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The field amplitude is considered here as a measure of 
the degree of excitation of a surface plasmon located on 
a specific interface. Figl^L shows the local amplitude of 
resonant orders (±1,0) at the substrate/metal interface. 
When the light crosses the substrate/metal interface at 
the first place, it shows a peak amplitude at the wave- 
length (3), predicted by EqQ] One notes a second peak 
at wavelength (3') close to (3), possibly related to inter- 
ferences of a true resonance with a non-resonant process, 
as found in Fano resonances^. These peaks can easily be 
associated with the dip at similar wavelengths in the re- 
flection spectrum. By contrast, when the vacuum/metal 
interface is illuminated first, the amplitude of these res- 
onances collapses with the result that there is an atten- 
uation of the associated dip in the reflection spectrum. 

FiglSb shows resonant orders (±1,0) on the vac- 
uum/metal interface. Here the resonance develops at 
the wavelength (1) and we note the existence of two 
other peaks at wavelengths (1') and (2'). It is likely that 
these peaks arise from a surface plasmon splitting due 
to the coupling of the (±1, 0) plasmon orders at the vac- 
uum/metal interface with the (±1,±1) plasmon orders 
at the substrate/metal interface. If the vacuum/metal 
interface is illuminated first, amplitudes of peaks (1) and 
(1') are stronger than in the opposite illumination. This 
explains the nonreciprocity found in the reflection spectra 
for wavelengths close to (1'). The spectral region close 
to the wavelength (2'), of the resonant order amplitude, 
does not significantly change when changing the propaga- 
tion direction, leading to a weak nonreciprocity effect in 
the reflection spectrum. In Fig. 3c we show the amplitude 
of orders (±1,±1) at the substrate/metal interface and 
focus on the resonance wavelength noted (2). A major 
contrast is obtained, with a large fleld amplitude under 
the direct illumination of the substrate/metal interface 
and a much weaker amplitude under illumination from 
the vacuum side. 

The amplitude of the resonances of orders (±1,±1) 
at the substrate/metal interface and (±1,0) at the vac- 
uum/metal interface are influencing each other because of 
the proximity of their resonance wavelengths. Indeed, the 
profile of the orders (±1,0) amplitude at vacuum/metal 
interface, just after the wavelength (1), can be under- 
stood as a Fano profile. Such a Fano resonance can be in- 
terpreted as resulting from the interference of the orders 
(±1, 0) resonance at the vacuum/metal interface, and the 
orders (±1,±1) resonance at the substrate/metal inter- 
face which then appears as a result of multiscattering of 
(±1,0) into (±1,±1). This explains, for instance, why 
the wavelengths (1') and (2') do not strictly correspond 
to the resonance wavelength of the different resonant or- 
ders. Fano profiles resulting from the interference of each 
order as a result of multiscattering lead to a visible res- 
onance wavelength shifts. 

We note that this analysis of the comparison between 
surface plasmon features and reflection corroborate the 
interpretation of Altewischer et a&, from which useful 
information about surface plasmons can be drawn from 



the study of the reflection nonreciprocity. 

In addition to transmission and reflection, it might be 
useful to investigate the loss response function for both 
directions of illumination. This is shown in Fig|31i. As 
light comes from the vacuum side, we observe two sharp 
peaks at wavelength (1') and (2'). These peaks corre- 
spond to the surface plasmon resonance already described 
above. On the other hand, when the light enters through 
the glass side, the first peak collapses, whereas another 
large double peak appears at the wavelength (3'), in 
perfect correlation with the surface plasmon excitation. 
This, again, is in full agreement with the experimen- 
tal results of Barnes et a'fi and previous computations^, 
confirming that both loss peaks correspond to a surface 
plasmon resonance. The loss line shapes can easily be 
correlated with the dominant features of the reflection 
spectra. These features indicate that, when SPs are ex- 
cited, a large part of their energy is dissipated via the 
Joule effect, leading to lower reflection intensities. In ad- 
dition, a smaller part of SPs energy is transferred to the 
transmitted beams as a result of multiscattering. So, as 
SPs are essentially excited on the interface which is en- 
countered flrst, the reflection spectrum is dominated by 
processes associated with the losses from the plasmons 
carried by this interface. 

As already noted by Barnes^, this is quite different for 
transmission, where SPs of both interfaces equally con- 
tribute to enhanced transmission whichever side is illu- 
minated first. The interface which receives the incident 
beam provides strong SPs resonances. Such resonances 
only give a weak contribution to the transmission since 
the associated field is damped when propagating through 
the holes before providing transmitted intensity. On the 
opposite interface, SPs are excited by an incident light 
weakened by its propagation through the holes, so that 
SPs on the opposite interface are only weakly excited. 
Their contribution to the transmitted field is thus also 
weak. In other words, SPs of the illuminated interface 
are strongly excited and then damped after propagation 
through the holes, whereas SPs on the opposite interface 
are weakly excited as the incident light is first damped. 
The result is that SPs of both interfaces contribute to 
the transmission with similar magnitudes and that the 
choice of the direction of propagation does not infiuence 
the information content of the measured transmission. 

In conclusion, the examination of several experimental 
studiesii in the light of the interpretation scheme pre- 
sented elsewhere^ reinforces the point of view expressed 
in this scheme by exhibiting qualitative and quantitative 
agreements between measured and computed light filter- 
ing spectra of a perforated gold layer covering a glass 
substrate. A particular effort was made to study the non- 
reciprocity property of the reflection. It was found that 
the asymmetry relative to the light propagation direction 
is effectively caused by the localization of SPs resonances 
at gold film surfaces and that this can be exploited to 
investigate separately the gold/glass and gold/ vacuum 
plasmon coupling to radiative waves. The transfer func- 
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tions computed in detail in this paper agree well with 
the quoted experimental data, available to us at the time 
of writing. One of the outcomes of these experiments is 
the role played by the corrugated-surface plasmons in the 
build up of the light transport process. In order to shed 
some light on this issue, the present work also looked 
into the theoretical value of the excited field at the loca- 
tion of the metal/ vacuum and metal/glass interfaces and 
determined the degree of excitation of the interface plas- 
mons. The strong correlation found between the plasmon 
amplitude and the reflection variations shows unambigu- 
ously the role played by the metal surface collective ex- 
citations in the light transfer process. This conclusion 
suggests that the experiments could be taken some steps 
further by considering complementary measurements, at 
least on the outer vacuum/gold surface, of the far field in 
reflectance experiments and the near field with the scan- 
ning near-field optical microscope. A correlation between 
these separate investigations should be highly instrumen- 
tal in elucidating the mechanisms of the light transfer 
processes. 
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